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■ ABSTRACT 

o 

We propose a new two-stage model for acceleration of electrons in solar flares. In 
D ■ the first stage, electrons are accelerated stochastically in a post-reconnection turbulent 

downflow. The second stage is the reprocessing of a subset of these electrons as they 
pass through a weakly compressive fast shock above the apex of the closed flare loop on 
their way to the chromosphere. We call this the "shock reprocessing" model. The model 
reproduces the energy dependent arrival time delays observed for both the pulsed and 
£T) ' smooth components of impulsive solar flare x-rays with physically reasonable parameters 

m 

Q\ ■ as seen in the Masuda flare. The loop-top source distinguishes the shock reprocessing 

o 
in 

pendence of footpoint pulse strengths and the location and spectrum of the loop-top 
emission, and can account for the observed soft-hard-soft trend in the spectral evolu- 
tion of footpoint emission. Our model highlights the concept that reconnection is an 
acceleration environment rather than a single process. Which combination of processes 
operate may depend on the initial conditions that determine, for example, whether the 
reconnection downflow is turbulent. The shock reprocessing model comprises one such 



combination. 



1. Introduction 

An important observational constraint on the acceleration and transport processes in solar 
flares is imposed by X-ray arrival time delay measurements (Aschwanden, et.al. 1995; Aschwanden 
et al. 1996a,b, 1997, 1998; Aschwanden 1998; Aschwanden et al. 1999). It is observed that the 
non-thermal hard X-ray emission consists of two separable components, a smooth, slowly varying 
background, and a pulsed, rapidly varying modulation. Detailed studies of the temporal structure 
of the two components show that they possess several opposite characteristics: Typically, the 
variations in the smooth component are first observed at lower energies, then at increasingly higher 
energy. The arrival time delay between 25keV and 150 keV electrons is of order 5s. Conversely, 
the pulses are observed first in the highest energies, with lower energy photons appearing later; the 
lag time for 25keV photons relative to 150keV photons is of order 50ms. The trap-precipitation 
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model of Melrose & Brown (1976) offers one possible explanation for the time delay measurements. 
In this model, electrons are accelerated to high energy above the loop in the turbulent downflow 
and are then injected into the closed flare loop on the pulse time scale. Some electrons are trapped 
between magnetic compressions in the loop before precipitating onto the chromosphere, while the 
remainder precipitate directly. The two populations produce the smooth and pulsed components 
respectively. 

Although the trap precipitation model may plausibly account for the time delay observations in 
some flares, some significant uncertainties remain. While the needed injection may arise naturally 
as a feature of many acceleration models which are situated in the reconnection current sheet 
(i.e. Litvinenko (2000)), it is not clear how injection would take place in models which employ 
acceleration in the turbulent downflow region which may form between the reconnection sheet and 
the closed soft x-ray loops. In addition, Larosa & Shore (1998) argue that the pulsed structure 
can be produced by fluctuations in spectral hardness instead of a pulsed injection of non-thermal 
electrons. Furthermore, the model predicts that coronal X-ray emission emanates from within the 
trap region. The trap precipitation model places the trap, and thus this emission site, within the 
closed flare loop itself. This is seemingly at odds with at least some observations, such as those 
of of Sui et al. (2004); Masuda, et.al. (1996); Alexander & Metcalf (1997) which reveal a compact 
non-thermal source above the loop-top during the main phase of some flares. Another disadvantage 
of the trap precipitation model is that it does not retain temporal information of the acceleration 
process; such information is lost at the injection point. A different scenario in which the properties 
of the acceleration mechanism are more closely linked to the emission may provide more opportunity 
to test both the time delay mechanism and the acceleration mechanism observationally. We propose 
a new model called the "shock reprocessing model" as just such an alternative; it circumvents the 
loop injection problem, places the coronal emission above the loop-top, and retains some details of 
the acceleration mechanism in the footpoint emission. 

The shock reprocessing model employs both second order stochastic Fermi acceleration (STFA) 
and first order acceleration at fast shocks. Electrons are accelerated via STFA in the turbulent 
downflow region located below the reconnection point. Due to the fairly long acceleration time, 
typically a few seconds, variations on the eddy time, ^ Is, are smoothed out in the temporal 
spectrum of the STFA accelerated population. STFA, operating on a time scale of a few seconds, 
produces the smooth component of the X-ray emission described in Aschwanden, et.al. (1995). 
As the flow impinges upon the flare loop, typical flow speeds are marginally super-magnetosonic 
(e.g. Blackman & Field 1994). A weakly compressive fast shock forms. Because the flow is only 
marginally super-magnetosonic, the shock can disappear and reform as the flow speed fluctuates on 
the eddy time. This is a key feature of the shock reprocessing model. First order Fermi acceleration 
occurs at the fast shock and the transient nature of the shock can be responsible for the pulsed 
emission structure observed in impulsive solar flares. Furthermore, the shock is strongest above the 
apex of the loop, and diminishes toward the wings. The shock formation condition is typically met 
over a portion of the loop; as a result, only a fraction of the STFA accelerated electrons undergo 
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the second phase of acceleration.. 

We discuss the details of the acceleration mechanisms in section 2. In section 3 we determine 
the range of shock parameters required by the shock reprocessing model, as well as a simple model 
of shock formation. In particular, we calculate the compression ratio and spatial filling fraction of 
the shocks. In section 4 we discuss the role of cooling in the STFA region in matching the observed 
smooth component time delays and the presence of loop-top emission sites in a small fraction of 
Yohkoh flares, such as the Masuda flare (Masuda et al. 1994). The distinct presence of these X-ray 
sources just above and exterior to the loop-top is a distinguishable feature of the shock reprocessing 
model when compared to the trap precipitation model. We conclude in section 5. 

2. Flare properties and acceleration mechanisms 

Before examining the shock reprocessing model in depth, we briefly review models of flare mor- 
phology and some properties of STFA and shock Fermi acceleration which will later be employed. 

2.1. Flare properties 

It is well understood that solar flares are driven by magnetic reconnection in the solar corona 
(i.e., Priest & Forbes (2002)). The reconnected field lines relax into a closed loop structure evi- 
denced by the soft x-ray loops which are filled with hot plasma. Hard x-rays are observed at the 
chromospheric footpoints of the closed loops as well coronal hard x-ray sources in some flare. The 
acceleration mechanisms responsible for the hard-xray producing electrons are not yet fully known; 
candidates include DC field acceleration (Litvinenko 2000), stochastic acceleration (Miller, Larosa, 
Sz Moore 1996; Chandran 2004; Selkowitz & Blackman 2004), and shock acceleration (Tsuneta &; 
Naito 1998), or combinations of these processes (Somov & Kosugi 1997). Further complicating the 
issue is that the downflow plasma between the reconnection point and the relaxed loops may either 
be laminar or turbulent. Laminar flow would favor stationary (Melrose k, Brown 1976)and col- 
lapsing Karlicky & Kosugi (2004) trap precipitation models of electron acceleration, which require 
an ordered field structure. Turbulent downflow would favor stochastic acceleration models, which 
require turbulent field structure. 

In addition, there are two possible directions of motion to consider for the evolution of the 
reconnection point. The standard model (i.e. Priest & Forbes (2002)) focuses on the vertical 
direction of motion in which the reconnection point moves upward through the corona over time. 
As a result, the footpoint separation and height of the soft x-ray loop increases with time. However, 
recent observations with RHESSI indicate that at some flares have a "zipper" morphology; the 
reconnection evolves laterally along an arcade structure at near constant height, with the soft x-ray 
loop following the reconnection point along the arcade at near constant footpoint separation and 
loop height (Grigis & Benz 2005). Regardless of the direction of motion of the reconnection point, a 
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2-D cross section of flare morphology is shown in figure 1. The region labelled outflow may be either 
turbulent or laminar, depending on the local Alfven speed and the initial gradient in the downflow 
across the downflow cross section (Chieuh & Zweibel 1987). The shock reprocessing model focuses 
primarily on the vertically evolution with turbulent downflows. 

It is important to note that there is also a variety of coronal hard x-ray emission characteristics 
observed. The discovery event, dubbed the Masuda flare (Masuda et al. 1994), shows hard x-rays 
above the closed loop with a non-thermal kernel (Alexander & Metcalf 1997) as well as a very hot 
thermal region. On the other hand, many coronal sources appear to be wholly thermal (Emslie 
2003). Some are located inside the soft x-ray loops, not above (Veronig & Brown 2004). The variety 
of x-ray emission properties and reconnection morphologies implies that even if all solar flares result 
from a basic reconnection environment and a soft x-ray loop, they are not explained within a single 
scenario of specific processes such as trap precipitation or shock reprocessing. Key elements, such 
as laminar vs. turbulent downflow may determine the differences between some classes of flares. 
A necessary condition for turbulence to form in the reconnection outflow, super- Alfvenic downflow 
speeds, is also an approximate result of reconnection outflows (Blackman & Field 1994; Karlicky 
& Kosugi 2004). Ultimately, models which consider flares with and without turbulent downflows 
need to be comparatively explored. 



STFA is capable of accelerating solar flare electrons out of a thermal population (Miller, Larosa, 
& Moore 1996; LaRosa et al. 1996; Chandran 2004; Selkowitz & Blackman 2004) and up to high 
energies with a spectral shape that depends mainly on the mechanism which traps electrons in the 
acceleration region. In the reconnection outflow, trapping is provided by wave particle interactions 
which also serve as the pitch angle scattering necessary to sustain acceleration (Achterberg 1981). 
Below a transition energy Et ^ lOkeV, pitch angle scattering in solar flares occurs through interac- 
tions with whistler waves. In the presence of these waves, STFA produces a quasi-thermal electron 
spectrum below E t . Above E t , the scattering mode required to produce a power law spectrum is 
constrained but undetermined. This remains a key unresolved issue for STFA. In order to produce 
a power law, the pitch angle scattering length must be inversely proportionl to the electron energy 
E (Selkowitz & Blackman 2004) 



where X p is the pitch angle scattering length scale, L is the linear size of the acceleration region, 
and T is a proportionality constant which is fixed by the scattering physics. The resulting spectrum 
has the form 



2.2. Power law acceleration solely by STFA 
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where N(E) is the number of electrons with kinetic energy E, and the spectral index j st is given 
by 

where m p is the proton mass and va is the local mean Alfven speed in the plasma. Variations in T, 
resulting from local changes in the pitch angle scattering length scale, produce the observed range 
of spectral indexes among flares. Typically, STFA electron (X-ray) spectra have a mean index of 
4.5(3.5) (Bromund, McTiernan & Kane 1995). 



2.3. Power law acceleration solely by fast shocks 

Shock Fermi acceleration had been considered as a means of producing the power law elec- 
trons required to explain flare X-rays, for example Blackman (1997); Tsuneta & Naito (1998). A 
stationary fast shock can form at the point where the turbulent reconnection outflow plows into 
the loop. This shock, located just above the top of the loop, can accelerate electrons into a power 
law distribution. 

Shock Fermi acceleration is well studied (Bell 1978; Bell 1978b; Blandford & Eichler 1987). 
In the standard theory, charged particles gain energy from repeated transits across the shock. 
Since there is a net energy gain in each cycle, the process is first order and results in a rapid 
acceleration. The energy spectrum produced by shock Fermi acceleration in impulsive solar flares 
in the non-relativistic regime is a power law of the form 



N(E)ocE- Sr , (4) 

where 5 r = 2(r + 2)/(r — 1), and r is the compression ratio across the shock (e.g. Jones & 
Ellison 1991). Furthermore, the shock compression ratio is related to the downflow Mach number 

M = Vf/c s 

( 7 + l + (7-l)(M 2 -l))' 1 ' 

where 7 is the adiabatic index of the plasma. From simple theoretical estimates, it is expected 
that the shocks are weakly compressive (Blackman & Field 1994; Blackman 1997). To fit the 
observationally inferred electron spectrum (Brown 1971; Bromund, McTiernan & Kane 1995) of 
(S r ) = 4.5 requires a compression ratio r = 2.6, significantly less than the maximum of r = 4 for a 
non-cooling shock, and consistent with theoretical predictions. 

Weak fast shocks are expected, but standard shock Fermi acceleration also requires injection 
of high energy electrons which satisfy two conditions: their gyro-radius must be larger than the 
shock thickness (roughly the thermal proton gyro-radius) to see the shock as thin, and the electrons 
must be energetic enough to scatter off of upstream Alfven waves (Bell 1978; Blandford & Eichler 
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1987). The threshold energy, Ei n for injection at impulsive solar flare shocks is found to depend on 
the angle of obliquity between the normal to the shock surface and magnetic field lines, and ranges 
from 2keV at 85deg to lOkeV at Odeg (Tsuneta & Naito 1998). Herein we propose that STFA, in 
the presence of pitch scattering by whistler waves, is the shock acceleration injection mechanism; 
it accelerates electrons from the background plasma to energies above the proton thermal energy. 
Note that this scenario differs from that of Tsuneta & Naito (1998) in which whistler waves are also 
called upon as a scattering mechanism, but the acceleration is single stage; all electron acceleration 
takes place at the shock in their model. 

The electron spectrum from fast shock acceleration following STFA in the presence of whistler 
waves is can be reproduced in some STFA models without shock acceleration, although the spectrum 
from shocks depend only on the compression ratio. However, in the case with shocks, the transition 
to a power law spectrum occurs at the shock injection energy, Ei n , instead of at E%. A similarity 
between STFA and shocks is that both have difficulty in producing downward spectral breaks from 
an initially steep spectrum but can produce upward breaks (hardening). It has been established 
by Bell (1978b) for example, that processing of an electron population by a shock can harden an 
already existing power law distribution, with spectral index ji n , or reduce the index below ji n , but 
cannot soften it, or raise the index above 7j n . If there is a power law spectrum downstream of a 
shock with a low energy cutoff below the injection energy and the shock is sufficiently strong, the 
upstream spectrum will be harder above Ei n but remain unchanged below Ei n . This property is 
important to the shock reprocessing model. 

3. The shock reprocessing model 

In the following section, we present the shock reprocessing model and compare it to the trap 
precipitation model of Brown and Melrose (1976). 

3.1. Overview of the model 

Observations of non-thermal X-ray emission at the footpoints of solar flares indicate that the 
observed emission can be divided into two components, a smoothly varying part and a rapidly 
varying set of pulses (Aschwanden, et.al. 1995). Observations in multiple energy bands show that 
for the smooth component, variations in X-ray intensity arrive first at lower energies, and later at 
higher energies. The pulses exhibit an opposite time delay structure; variations are seen first in 
higher energy bands. The series of papers by Aschwanden and collaborators (Aschwanden, et.al. 
1995; Aschwanden et al. 1996a,b, 1997, 1998; Aschwanden 1998; Aschwanden et al. 1999) models 
the X-ray arrival time observations in solar flares within the trap precipitation model. We propose 
the shock reprocessing model as an alternative. 

In the trap precipitation model (Fig. 2), reconnection occurs high above the flare loop. Elec- 
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trons are then accelerated to high energies above the loop over a short time scale, resulting in 
a pool of non-thermal electrons in a near equilibrium distribution. An injection event somehow 
loads a pulse of electrons, with an energy spectrum matching that of the pool above the loop-top, 
onto the loop. Injection occurs over a very short time, 50ms, matching the observed X-ray pulse 
durations. Injected electrons travel down the flare loop, where they encounter a magnetic mirror 
as the loop field lines converge. Electrons with sufficiently high pitch angle relative to the local 
field are trapped at the mirror, while those with small pitch angle stream through. This divides 
the electrons into two populations: directly precipitating, and trap precipitating electrons. The 
directly precipitating electrons do not reflect at the mirror, and proceed directly to the dense foot- 
point region. There they emit X-rays via thick target bremsstrahlung; these electrons retain the 
short pulse time scale of the injection, and thus produce the pulsed emission component. Since 
the pulse electrons are injected simultaneously at all energies, the main time dispersion of a pulse 
results from time of flight down the loop, with higher energy (faster) electrons arriving earlier than 
lower energy (slower) electrons, consistent with the pulse observations. The second population, 
trap precipitated electrons, are mirrored at the magnetic field compressions due to their large pitch 
angles. These electrons remain trapped between the mirrors on opposite sides of the flare loop 
until such time as their pitch angles scatter to sufficiently small values that they can pass through 
the mirrors. Aschwanden et al. (1997) find that the slow time delays are consistent with loop-top 
trapping, where escape from the trap is governed by the collisional scattering time which varies 
as E 3 / 2 . Because the trapping time is significantly longer than the injection time, electrons from 
many injections escape together, creating a smooth component. As a result of the positive energy 
dependence of the trapping time, lower energy electrons escape the trap earlier than higher energy 
electrons, and the smooth emission component exhibits variations at lower energies first. The trap 
precipitation model offers an explanation for the two component emission, as well as the time delay 
data. However, it is not unique. 

The shock reprocessing model is an alternative to trap precipitation. An outline of the model 
is shown in Fig. 3. As in trap precipitation, the shock reprocessing model begins with reconnection 
high above the flare loop, followed by acceleration above the loop-top. However, in the shock 
reprocessing model the initial acceleration is slow and consistent with STFA in the reconnection 
outflow, occurring on a time scale of 5s (Chandran 2004; Selkowitz & Blackman 2004). Following 
STFA, the electrons flow into the region just above the loop-top. At or just above the loop-top a 
weakly compressive fast shock forms (Blackman & Field 1994; Tsuneta 1996). The shock can vary 
in spatial extent and compression ratio. A fraction of the STFA accelerated electrons also pass 
through the shock and are accelerated a second time (reprocessed) . The remainder pass directly to 
the flare footpoints. 

Like the trap precipitation model, the shock reprocessing model involves two populations of 
electrons, but in this case they are the directly streaming (STFA only) electrons and the shock 
reprocessed electrons. The directly streaming electrons produce the smooth X-ray component. 
Since there is no trapping or injection within this model, variability in the produced emission 
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results solely from the acceleration process. The ~ 5 s acceleration time for STFA (Selkowitz & 
Blackman 2004; Chandran 2004) is consistent with the delay time observed between 25keV and 
150keV electrons in flare observations (Aschwanden, et.al. 1995; Aschwanden et al. 1996a, b, 1997, 
1998; Aschwanden 1998; Aschwanden et al. 1999). The reprocessed electrons are those which pass 
through the shock and undergo a second, rapid acceleration. Subsequently, they travel without 
coronal trapping to the footpoints. Because the shock acceleration time scale is extremely short, 
reprocessed electrons can be considered effectively injected at the loop-top simultaneously across all 
energies. The number of electrons which are reprocessed is determined by the spatial extent of the 
shock. We proceed by developing a phenomenological procedure for determining pulse strengths, 
discussing shock formation and geometry as further constraints on the model, and predicting the 
location of the coronal emission site coincident with the STFA and cooling in a region above the 
loop-top. 

The absence of an above the loop-top, or loop-top trap is a key distinction between shock 
reprocessing and trap precipitation models. Expansion of the closed loop is often assumed to occur 
near the loop-top and is a requirement of the trap precipitation model to create the trap in the 
form of a magnetic mirror. Recently however, Bellan (2003) presents observational and theoretical 
evidence against the formation of a trap, finding instead that the loops are very nearly axially 
symmetric from footpoint to apex. This poses a challenge to the static trap models which require 
narrowing loop cross sections at the footpoints and broader loop cross sections at the apex. 



The filling fraction of the shock, F, can be understood as both the fraction of the downflow 
cross section which encounters the shock, and as the fraction of non-thermal electrons which are 
reprocessed. Note that this sets the firm upper bound F = 1. The physics relating F to shock 
formation is deferred to section 3.4. Here, we constrain F from observations and power law models 
of the electron and x-ray spectra. We begin by assuming that the electron spectrum resulting from 
STFA is a single power law, given by 



where N s is the total number of electrons per unit energy per unit time passing through a slab of 
unit surface area coplanar with the shock (henceforth we refer to similar quantities as areal fluxes), 
and is therefore the total non-thermal electron population, as a function of energy. Nq s is the areal 
flux of electrons at energy Eq s , the lower threshold energy for the onset of the STFA power law 
spectrum, and 5 S is the spectral index of the STFA electrons. The total number of non-thermal 
electrons is then given by 



3.2. Parameterizing the model 




(6) 
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= aT N sd E = ^^, (7) 

J En, O s -l 



where Nt is the total number of non-thermal electrons and A is the area of the downflow. Because 
of the injection criterion, shock reprocessing does not increase the total number of non-thermal 
electrons. As a result, the total number of reprocessed non-thermal electrons, N p , can be written 
as 



— fa r N , dE=F ^m- ls -- 1] =FE j!i>)-' s --\ T , (8 ) 

JE r d s -l \E 0s J \ E 0s J 

where we have taken into account only those electrons with energy above the shock injection 
threshold, Eo r , and assume Eq t > Eq s . 

We can obtain a second expression for N p by assuming a power law of similar form to Eq. (6) 
for the reprocessed electrons, 

(E \ 

where quantities with the subscript r are defined similarly to those with an s in Eq. (6), but for 
the shock reprocessed population. The total number of electrons in the reprocessed population, 
N p , is obtained analogously to Eq. (7) by integrating over energy 

N p = FA r N r dE = p^E^L. (10) 

JEor r - 1 

Setting Eqs. (8) and (10) equal yields 

For thick target bremsstrahlung at the footpoints, the photon spectrum is likewise a power law 

(E \ ~ ax 
IF- ' ( 12 ) 

where x can be either r or s for the shock and STFA components, M x is the areal photon flux at 
energy E, and Mq x is the areal photon flux at energy Eq x . The spectral index for the photons is 
related to that of the electrons by a x = 5 X — 1 (Brown 1971). Combining this with Eqs. (6) and 
(9) yields the additional constraint 

^- 1 >tr (< - 1 >lr (13> 
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Integrating both the STFA and reprocessed spectra of Eq. (12) from the minimum observed energy 
Emin, where we assume E m i n > Eq t > Eq s , gives the areal photon flux for both the STFA and 
reprocessed components. 

M t = aT M s dE = A^f (^-Y iSS ~ 2 \ (14) 



MftrEftr I Er, 



-{Sr- 2) 



M p = Al Mr dE = A^f[^\ , (15) 

JE min d r -2 \ E 0r J 

where Mp and M p are the energy integrated areal photon flux from the STFA and shock reprocessed 
populations respectively. Using Eqs. (11), (13), (14), and (15) we can rewrite M p in terms of Mp 



M v=[ 5 -f-^)El- 5 °E^ Ml , (16) 



Given the above, we can construct photon fluxes for the pulsed and smooth X-ray components 
given values of the free parameters of the model: F, 5 r , and 5 S . We define the smooth and pulsed 
emission as follows. The smooth emission is the total rate of photon counts at the footpoints in 
the absence of pulses, which is simply given by Eq (14) times the area of the downflow. The pulse 
emission, Mp, is the enhancement above the smooth emission, which we take at the peak of the 
pulse, and is given by 

M P = FM P -FM T = F ((|^|) E 5 ;- S °E 6 ^ _ ^ Mr> (17) 

which is the total emission by the reprocessed electrons, FM p , less the emission that would have 
been produced by those electrons if they were not reprocessed, FMp- For further simplicity, we 
choose to measure energy in units of Eo r . Thus, we can write 

u-£-f((^)ei!s+-i). (18) 



Mp \\S r -2 / 

The model produces pulses of magnitude M = Mp/Mp determined by four parameters: the shock 
filling fraction F, the STFA electron spectral index <5 S , the reprocessed electron spectral index 5 r , 
and the minimum observed energy E m i n . 



3.3. Fitting the observed pulse strengths 

Fig. 4 shows a set of curves of constant pulse strength (M = 0.1,0.2,0.4,0.6) in 5 r vs F 
space. We have taken E m i n = 5Eo r /3 and S s = 4, consistent with the low energy observational 
cutoff of E = 25keV in Aschwanden, et.al. (1995) and a presumed non-thermal cutoff energy of 
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^ 15keV. A number of important properties of the model are evident. First, the curves diverge 
asymptotically toward F = oo as S r approaches 5 S . This is expected from both Eq. (18) and 
the physics of shock acceleration. The result of Bell (1978b) that a shock can harden, but not 
soften, an electron spectrum indicates that the shock reprocessing model can only produce pulses 
via reacceleration at the shock if S r < 5 S . Since 5 r is determined by the shock compression ratio, 
pulsing (as opposed to steady enhancement) of the reprocessed component requires variation of 
the shock compression ratio, or equivalently the downflow Mach number, on the pulse time scale. 
A second property of the curves is that there is an upper limit on 5 r for any given observed M, 
which is set by the F = 1 line. By definition, F > 1 represents an unphysical solution, where 
the number of reprocessed electrons exceeds the total number of non-thermal electrons. The shock 
reprocessing model thus can constrain parameter space. Any choice of M and S r which predicts 
F > 1 represents a physically unrealizable state. For example, Fig. 4 shows that pulses stronger 
than M = 0.6 cannot be formed by shocks with 5 r > 3.4, while pulses of strength M = 0.4 
can be formed out to 5 r = 3.6. Generally, the stronger the pulse, the harder the limiting shock 
spectrum. Finally, notice that all of the curves converge toward F = at 5 r = 2. This results from 
the increasing difference between the reprocessed and STFA spectra as the shock approaches the 
limiting case of 8 r = 2 at a strong shock. At all spectral indexes, there is a high end cutoff photon 
energy set by the total energy budget of the flare. However, for all values of S r > 2, the integrals of 
Eqs. (14) and (15) are dominated by the lower energy bound, so ignoring the cutoff is acceptable. 

Fig. 5 shows the effects of varying the value of E m i n on the curves of constant M in F and 
8 r space. We have taken 5 S = 4 for each of these plots. As we raise E m i n , the curves shift toward 
smaller values of F for any given S r . The shift results from moving the observational cutoff energy 
further from the onset of the reprocessing power law. The difference between the two power law 
spectra, M r and M s is greater at higher energies. Moving the lowest energy of observation further 
away from the shock injection energy results in stronger pulses. Therefore, the model predicts 
that observed pulses will be stronger in higher energy bands for any given flare. We explore the 
ramifications of this prediction in section 4. Likewise, Fig. 6 shows the effects of varying 5 S while 
keeping £ m „ fixed. In this set of plots, E m i n = 5Eo r /3, and S s takes the values 3.5,4, 4.5, and 5. 
The asymptotic divergence of F always occurs at 5 r = 5 S . Increasing 5 S results in the formation of 
stronger pulses at any given values of 5 r and F. 



3.4. Constraints of shock formation 

An additional constraint on the parameter space available for shock reprocessing can be ob- 
tained from studying shock formation. A detailed study of this problem is beyond the scope of the 
current work, but a simplified treatment can provide some insight. Fig. 7 shows a schematic model 
of the shock forming region. To make the problem tractable, we assume the flow forms a shock at 
any point along the loop-top where the component of the flow speed normal to the loop surface is 
super-fast-magnetosonic Furthermore, we take the shock to be planar, and compute an averaged 



compression ratio along the shock. The loop geometry is taken to be elliptical, consistent with the 
measurements of Aschwanden et al. (1996a). The tangent line to the loop surface has a slope 



s = 




(19) 



where s is the slope of the tangent, h is the loop height, 2w is the footpoint separation, and x is 
the distance from the loop center. The component of the flow normal to the loop v± is given by 



where Vf is the downflow speed. Setting Eq. (20) equal to the fast-magnetosonic speed, cy, 
substituting for s from Eq. (19), and solving for x gives the critical distance, x max , for shock 
formation. The covering fraction of the shock is given by F = x max /w. The compression ratio of 
the shock is given by evaluating Eq. (5) at each point in the range < x < x max and averaging. The 
resulting spectral index 5 r of electrons accelerated at the shock obtains from Eq. (4). Solutions of 
these equations are plotted in Fig. 8 for various values of h/w. For all loop heights, F approaches 
1 as the flow approaches the strong shock limit of S r = 2. As S r grows large, the flow velocity 
decreases, and F = at the limit of v = c/, or 6 r = oo. Notice that as the loop gets more elongated 
(larger h), F decreases. 

We overlay Fig. 8 onto panel (c) of Fig. 5 to produce Fig. 9, which incorporates the two sets 
of constraints: the pulse strengths from the shock reprocessing model, and the shock formation 
physics. The physically reasonable region, F < 1 and S r < 4, is delineated by the box. Selecting 
values for the two direct observables, loop height h and pulse strength M picks out a pair of curves, 
one from each set. These curves have a single intersection, fixing values of F and 5 r . Observations 
indicate that typically 0.1 < M < 0.6 and h/w ^ 1 (Aschwanden et al. 1996a). This places 
3 < S r < 5, consistent with the downstream fast shock spectral indexes calculated in Blackman Sz 
Field (1994) and the observed X-ray spectral indexes. 



An above the loop-top hard X-ray source has been observed in some impulsive phase flares, 
most notably the Masuda flare (Masuda, et.al. 1996; Tsuneta 1996; Sui et al. 2004), but not in 
others. During the main phase of the Masuda flare, a kernel within the thermal source is clearly 
observed to be non-thermal (Alexander & Metcalf 1997). It has been argued (Blackman 1997; 
Tsuneta &; Naito 1998) that the loop-top source is associated with trapping and acceleration of 
electrons below the stationary fast shock. In this scenario, the conditions for shock formation 
are met in a small fraction of flares, and then only marginally, resulting in weakly compressive 
downstream shocks in some flares, and no downstream shocks in others (Tsuneta 1996). The 
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3.5. Loop-top emission and smooth component time delays 
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low shock formation rate may explain the infrequent detection of non-thermal loop-top emission. 
It has also been argued (Petrosian et al. 2002) that the non-thermal loop-top source is actually 
prevalent in impulsive flares, but usually very dim, with typical intensity less than one tenth that 
of the footpoint sources. The contrast limit of Yohkoh is roughly 10; observations are limited 
to unusually bright loop-top sources or limb flares with obscured footpoint emission. We show 
that the shock reprocessing model is consistent with the latter scenario if STFA takes place in the 
above the loop-top emission region, and is accompanied by in situ bremsstrahlung cooling of the 
electron population. That being said, there are recent RHESSI observations of two flares which are 
dominated by coronal emission that is cospatial with the soft x-ray loop (Veronig & Brown 2004). 
We consider this latter issue further in section 5. 

Although the Masuda flare exhibits pulses, they are longer in duration than the footpoint 
pulses. Loop-top pulses are typically of duration > 10s (Mrozek &; Tomczak 2004), as opposed 
to the sub-second pulse structure in the footpoint emission. The pulsed component to which we 
attribute shock reprocessing is the sub-second pulse structure in the footpoint emission, not the 
much longer pulses which may come from a globally bursty reconnection. When the sub-second 
pulses are subtracted from the emission profile the remaining component we refer to as the smooth 
profile. In this respect, in addition to explaining the sporadic appearance of above the loop-top 
non-thermal sources, cooled STFA also produces the proper time delays for the smooth emission 
component. Aschwanden et al. (1997) analyzes the arrival time delays for the smoothly varying 
X-ray component of 78 flares. From these, they obtain plasma densities in the soft X-ray loops for 
44 events, rejecting the remaining 34. Of the rejected events, 29 were unsuitable to the analysis 
technique because the fast varying component could not be sufficiently deconvolved from the smooth 
component. The remaining 5 rejections were due to poor convergence of their fitting model. The 
remaining 44 flares have smooth component time delays consistent with the trap precipitation 
model, matching the predicted trapping time t(E) oc E 3 ^ 2 . 

This effect is also predicted by our shock reprocessing model. In shock reprocessing, the pre- 
shock STFA acceleration rate is expected to determine the smooth variation time delays. Variations 
in the smooth component are the result of changes in the STFA acceleration region, and appear at 
high energy on a time scale equal to the STFA acceleration time to that energy. The time delay 
curve should thus be obtained from the STFA acceleration rate: 



For non-relativistic STFA Q(E) is energy independent (Selkowitz & Blackman 2004) and r oc E. 
For statistical acceleration in the quasi-linear regime, Chandran (2004) finds the standard result: 



(21) 



where Q(E) is prescribed by STFA. We define 




(22) 
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Q(E) oc E, and r oc e E ^ E °. Neither of these forms is consistent with a good fit to r oc E 3 / 2 . Cooling 
in the acceleration region, which is neglected by Selkowitz &: Blackman (2004) can be invoked to 
reduce the acceleration rate. 

Consider the case where STFA with in situ cooling by bremsstrahlung occurs in a region above 
the loop-top shock. The power loss due to bremsstrahlung for a single electron is given by (Rybicki 
& Lightman 1979) 

P B = l.lmoE^keVs- 1 , (23) 

where the dimensionless electron density nio = n/(10 10 cm -3 ), E\ = E/lkeV, and we have assumed 
a constant Gaunt factor Gff = 1. Notice the linear density dependence for the Bremsstrahlung 
rate, which arises because we consider the mean cooling rate for a single electron in a background of 
protons at density n 10 (not the emission from an ensemble of electrons in a background of protons). 
Q(E) is thus given by 

Q(E) = (^f) -P B = W0n% 2 - l.ln 10 El /2 , (24) 
\at J stfa 

where we have adopted the STFA acceleration rate of Selkowitz & Blackman (2004), which reex- 
amined the STFA process and determined that, due to averaging over the limited range of pitch 
angles which are subject to reflection, as opposed to all pitch angles, the mean acceleration term 
dominates over the diffusive random walk component in the sub-relativistic regime. Taking the 
flare electron density to be nio = 1 an d performing the integral in equation 22 yields the function 
shown in Fig. 10 where we plot t{E) for both cooled and uncooled STFA as well as the empirical fit 
of r oc E 3 / 2 . Cooled STFA is consistent with the time delay measurements of Aschwanden, et.al. 
(1995); Aschwanden et al. (1996a,b, 1997, 1998); Aschwanden (1998); Aschwanden et al. (1999) 
over the observed energy range of 10 — 150keV, provided we ignore collisional energy transfer to the 
protons. We follow Aschwanden (1998) in taking such collisional losses to be negligible at energies 
above ^ 30keV, where the energy loss time is long compared to the travel time of an electron from 
the acceleration region to the footpoints. We now discuss Coulomb collisions further. 

It has previously been noted (Brown 1971; Schatzman 1965) that Coulomb losses due to 
electron-electron scattering are typically greater than the Bremsstrahlung cooling rate at hard 
x-ray energies in flare plasmas. While this is true, within the STFA acceleration region, the elec- 
trons spectrum is nearly thermal. The non-thermal post-acceleration spectrum results from the 
escape rate of electrons from the acceleration region (Selkowitz & Blackman 2004) and thus the 
escaped electrons. For the distribution within the acceleration region, Coulomb scattering among 
electrons does not strongly affect a given electron's acceleration time; on average, each electron 
receives as much as much energy due to Coulomb scattering as it loses. Coulomb scattering with 
protons is, however, a potential cooling mechanism, since solar flare protons are sub-Alfvenic and do 
not participate in STFA. The Coulomb loss rate for electron-proton scattering is given by (Spitzer 
1956) 
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= 2.5 x 10 -2 ni £? 1 1/2 ■ (25) 

ep 

The electron-proton loss rate is significantly lower than the Bremmstrahlung loss rate, and electron 
cooling within the limited STFA acceleration region is dominated by the latter process. This result 
holds because the plasma flow participating in STFA contains an already nearly thermal electron 
population which is not interacting with cooler electrons. In the surrounding area, which is filled 
with ambient plasma, STFA accelerated electrons which have escaped the acceleration region do 
interact with the surrounding cooler electrons, resulting in cooling of the electrons, reshaping of the 
spectrum, and heating of the plasma. This heated plasma may emit the thermal loop-top X-ray 
component observed in a number of flares. 

Energy loss due to loop-top Bremmstrahlung cooling is significant for high energy electrons 
where the total energy lost to cooling is greater than the electron energy when it leaves the accel- 
eration region. The energy lost to cooling is released in x-ray luminosity. This is a plausible source 
of the weak non-thermal loop-top emission regions observed in Yohkoh studies of flares (Masuda, 
et.al. 1996; Tsuneta 1996; Petrosian et al. 2002). 




4. Discussion of observational implications and predictions 

Three major observational constraints on any model of electron acceleration in impulsive solar 
flares are the production of proper pulse strengths, the energy dependent arrival time delays, and the 
appearance of loop-top emission sites in only a fraction of observed flares. The shock reprocessing 
model successfully accounts for all three of these features through a two stage acceleration process. 
In the first phase, electrons are accelerated via STFA with appreciable thin target bremsstrahlung 
cooling in the turbulent downflow region above the flare loop. The cooling of trapped electrons 
produces the loop-top source, which is often too dim relative to the footpoints to be observed by 
Yohkoh. Subsequently, a portion of the electron population undergoes diffusive shock acceleration 
at a weakly compressive stationary fast shock, splitting the electrons into two populations: the 
shock reprocessed pulse population, and the unshocked STFA population. Sub-second pulses result 
from modulation of the shock strength on short time scales. Since the predicted shocks are only 
weakly compressive, small changes in the reconnection outflow can remove or reestablish the shock. 

The presence of the shock below the loop-top emission site implies that the pulse structure 
observed in footpoint X-rays, and generated by variations in the shock compression ratio, are not 
present in loop-top sources. This prediction of the shock reprocessing model needs to be tested, and 
should be feasible over time as the RHESSI dataset grows to include large numbers of flares with 
detectable loop-top sources. Sui et al. (2004) find evidence of a smoothly varying loop-top emission 
source in three RHESSI observed flares. There is no apparent pulse structure in these sources. 
However, Sui et al. (2004) argue that these coronal sources are distinct from Masuda flare type 
loop-top emission sources. The observations of Mrozek & Tomczak (2004) show a two component 
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emission structure, with both pulsed and smooth emission, superficially similar to the footpoints. 
However, in these sources, the pulse time scale is typically 10s or larger, not sub-second as in the 
footpoints. The coronal pulses are attributed to modulations in either the reconnection or downflow 
environment. They are not inconsistent with the shock reprocessing model. Still, larger sampling 
statistics on variability in loop-top emission sources would be of great interest in determining the 
applicability of the shock reprocessing model to impulsive solar flares. 

An additional concern for any model of solar flares is the electron supply problem. For example, 
the hard x-rays emitted in the Masuda flare require a non-thermal electron throughput of 2 x 10 35 s _1 
(Masuda et al. 1994). The reconnection downflow for the same flare contained ^ 5 x 10 35 electrons 
s _1 (Tsuneta & Naito 1998). Thus the electron acceleration mechanism must either be highly 
efficient, accelerating a fraction of the downflow electrons of order unity, or a secondary supply of 
electrons must be available. STFA falls into the first category; the process is highly efficient, so 
the problem of supplying a sufficent flux of electrons is alleviated. Since the bulk plasma flows 
downward to the footpoints, including both the electron and proton populations, charge neutrality 
is maintained without the need to rely upon a return current. 

Shock reprocessing also makes predictions regarding the pulsed emission. Because the pulse 
component of the emission is caused by the harder shock reprocessed spectrum, the observed pulse 
magnitude is dependent on the energy bin, where higher energy bins show larger pulse strengths. 
From Eq. (18), M cx E^~^ r , where M is the relative pulse strength, and J? m j n is the minimum 
energy of the observations. Notice that raising E m i„ increases the pulse strength. In its simplest 
discussed herein, shock reprocessing predicts that observations in a sequence of energy bins would 
show stronger pulses at higher energies. Obtaining such data is within the capability of current 
instrumentation. The magnitude of the effect can be calculated by evaluating Eq. 18 multiple 
times, with E m i n taken to be the lower energy limit of each detection band. This provides the pulse 
strength for all x-ray emission at photon energy E > E m { n . To carry out the observational test, 
one canreformulate the definition of M slightly by reevaluating the integrals in Eqs. (14) and (15) 
with the upper limit E max set to the high energy end of the bin. Alternatively one may collectively 
bin all photons with energy E > E m i n for successively increasing values of E m i n . For simplicity, 
we perform the latter procedure, assuming energy bins with E m i n taken in 20keV increments from 
lOkeV up to 150keV. For all bins, E max = oo. The results are shown in table 1. The range of pulse 
strengths predicted for energy bins in observations of a single flare is very strongly dependent on 
the difference 8 S — 5 r . Even in the case of S r = 3.8, a large range of pulse strengths is expected 
0.07 < M < 0.54. This effect can be measured, and can be used to test the shock reprocessing 
model and constrain 5 r , even for moderately small values of F. 

Furthermore, the strong pulses at high energies are consistent with the soft-hard-soft (SHS) 
spectral pattern observed in many flare spectra, even down to subsecond time scales (Grigis & 
Benz 2004, 2005). Typically, the SHS pattern traces total luminosity. As flare emission cycles 
through pulses, the spectrum starts out soft at the low luminosity onset of the pulse, steadily 
increases up to a maximal hardness at the pulse peak, then decays back to a softer spectrum as 
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M 


E min (keV) 


S r = 3 


5 r = 3.4 


S r = 3.8 


10 


0.6 


0.26 


0.07 


30 


3.0 


1.1 


0.23 


50 


5.4 


1.7 


0.32 


70 


7.8 


2.2 


0.38 


90 


10 


2.6 


0.43 


110 


13 


3.0 


0.48 


130 


15 


3.4 


0.51 


150 


17 


3.7 


0.54 



Table 1: Pulse strengths M for increasing energy bins and a range of 5 r . 

the luminosity returns to the pre-pulse minimum. The shock reprocessing model generates pulses 
via strengthening of the shock compression, which corresponds to hardening of the reprocessed 
spectrum, as illustrated in Table 1. Small changes in the super-magnetosonic Mach number of 
weakly compressive shocks can produce significant changes in the compression ratio. For pulse 
strengths in excess of 1, the footpoint X-ray spectrum will have index 7 r , the shock spectral index, 
since the emission will be dominated by the reprocessed population. As the pulse progresses toward 
the peak luminosity, the shock compression ratio peaks as well, and so does the spectral hardness. 
During the decay phase, the shock compression decreases, and the spectrum softens, resulting in a 
SHS pattern. Moderate pulses, especially those with large F will still produce a weak SHS pattern 
simply due to the large number of reprocessed electrons. 

5. Conclusion 

Reconnection can be thought of as an environment for particle acceleration. Within the recon- 
nection region and ensuing outflows, a combination of acceleration mechanisms may operate. We 
have proposed the shock reprocessing model as one scenario of interest for solar flares when the 
downflows form reconnection sites are turbulent. 

Four major features of the shock reprocessing model developed in this paper are as follows: 

(1) The model posits a two stage acceleration process: STFA in a turbulent reconnection downflow 
followed by first order acceleration at a loop-top fast shock. Hard X-ray emission sources are 
predicted in the STFA region, as well as at the foot-points. 

(2) Pulsed emission is produced by variations in the compression ratio of the fast shock, and thus in 
the post-shock electron spectrum. The shock does not fill the entire cross section of the downflow, 
and is weakly compressive. It is seen that relative pulse strengths produced within the shock 
reprocessing model are reasonable for realistic values of the local plasma parameters, and a simple 
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model of shock formation. Furthermore, the generation of pulses at the shock predicts that the 
above the loop-top source does not exhibit the subsecond pulse structure observed at footpoints. 

(3) Both the fast (pulse) and slow (smooth) time delay measurements (Aschwanden, et.al. 1995; 
Aschwanden et al. 1996a, b, 1998; Aschwanden 1998; Aschwanden et al. 1999) can be explained 
within the shock reprocessing model. The fast time delays in the pulsed component are time of 
flight dispersion of the shock accelerated electrons. The slow time delays in the smooth component 
reflect the STFA acceleration rate. 

(4) The shock reprocessing model predicts an increase in relative pulse strength at higher energies. 
This prediction remains to be tested, and can distinguish between the shock reprocessing and trap 
precipitation models. 

The shock reprocessing model is not meant to be an exclusive solution to electron acceleration 
in all solar flares. It seems unlikely that the wide variety of flare phenomena can be explained if 
all reconnection regions have the same relative contributions for the combination of acceleration 
processes that can ensue. 

Generic features such magnetic reconnection, soft x-ray loops, and footpoint hard x-rays, do 
appear to be common to the vast majority of flares. However, other phenomena are not: Coronal 
x-ray emission sources appear sometimes with thermal and other times with non-thermal spectral 
characteristics, sometimes above the soft x-ray loop, or sometimes within it, and other times not at 
all. The reconnection morphology, while apparantly similar in any 2-D snapshot along the plane of 
the soft x-ray loops, can have either vertical (standard) or lateral ( "zipper" ) dynamical evolution. 

One particularly difficult observation to explain within the shock reprocessing model is the 
appearance, in a small number of flares, of coronal hard x-ray emission emanating from within the 
closed soft x-ray loop (Veronig & Brown 2004). These sources are non-thermal, and hard x-ray 
emission in the flares is dominated by the coronal, not footpoint sources. One might suspect such a 
source to be produced within the trap-precipitation model in the presence of an abnormally dense 
loop; the trapped electron population within the loop would emit via thick target Bremsstrahlung. 
Conversely, the shock reprocessing model does not load the electron beam onto the closed flare loop, 
and thus cannot produce these sources. This class of loop-top emission has only been observed in a 
few events. As the total number of coronal sources within loops grows, it would be of great interest 
to determine which, if any, environemntal parameters such as laminar vs. turbulent downflows 
correlate with the location of the coronal source. 

The phenomenological array of flares can likely be understood by a relatively small number of 
acceleration scenarios which each operate within the basic paradigm of magnetic reconnection and 
outflow. The power source driving all flares is reconnection high in the corona, which launches a 
downflow. Within the outflow, particle acceleration occurs, followed by x-ray, gamma ray, and radio 
emission in the lower corona and chromosphere. The particular acceleration scenario is determined 
by local environmental parameters in the reconnection, downflow, and emission regions. Shock 
reprocessing and trap precipitation are two such scenarios within this framework. A third scenario 
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seems to required to explain proton dominated flares (Hurford et.al. 2003; Miller and Roberts 1995; 
Miller, Emslie, and Brown 2004). 

In summary: the shock reprocessing model, whereby a fraction of stochastically accelerated 
electrons also passes through a weakly compressive stationary shock as they stream toward chro- 
mospheric footpoints, provides a scenario for explaining a variety of features in impulsive solar 
flares within the more general reconnection and outflow framework. The model produces pulsed 
and smooth spectral X-ray emission components consistent with the time delay observations of 
Aschwanden, et.al. (1995); Aschwanden et al. (1996a,b, 1997, 1998); Aschwanden (1998); Aschwan- 
den et al. (1999). Shock reprocessing also predicts the appearance of above the loop-top coronal 
emission sites observed in a fraction of flares. Furthermore, the model makes a sequence of testable 
predictions regarding the pulse strength as a function of energy and the SHS emission pattern. We 
have discussed how it is observationally feasible to test the shock reprocessing model. 
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Fig. 1. — The basic model of solar flare structure Blackman (1997). Reproduced by permission of 
the A AS. 
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The Trap Precipitation Model 
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Fig. 2. — A flowchart of the trap precipitation model, a is the electron pitch angle cosine and «o 
is the critical pitch angle cosine for trapping in the loop. 



The Shock Reprocessing Model 



Magnetic 
Reconnection 




Turbulent Downflow 
STFA on timescale 
t~5s 



Loop top 
Splitting into two 
electron populations 



passes through 
shock 



does not pass 
through shock 



Shock 

Reprocessed 
Electrons 



Directly 

Streaming 

Electrons 



Pulsed 
Emission 
at Footpoints 



7~ 



Smooth 
Emission 
at Footpoints 



Thick Target Bremstrahlung 



g. 3. — A flowchart of the shock reprocessing model. 
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(c) (d) 

Fig. 5. — Curves of constant M in S r — F space as E m i n is varied; 5 S = 4. (a) J? m j n = 3/3. (b) 
E m in = 4/3, (c) E m i n = 5/3. (d) E m i n = 6/3. E m i n is normalized to 15keV. In each panel, curves 
from left to right are M = 0.6, 0.4, 0.2, 0.1. Only F < 1 corresponds to physically realizable states. 
Increasing E m i n results in a lower required F for a given M 
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(c) (d) 

Fig. 6. — Curves of constant M in 8 S — F space as 5 S is varied, (a) 5 S = 3. (b) 5 S = 4, (c) 
S s = 5. (d) 5 S = 6. E m i n = 25keV in all four panels. In each panel, curves from left to right are 
M = 0.6,0.4,0.2,0.1. Only F < 1 corresponds to physically realizable states. Increasing 5 S implies 
a higher F is required to produce a given M. 
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Downflow 




w 

Fig. 7. — A schematic of the shock formation model. Vf is the downflow speed, v n is the component 
of the speed normal to the loop, s is the slope of the loop, w and h the loop width and height, and 
x the distance of the point at which s is evaluated from the mid-line of the loop. A fast shock forms 
above the loop-top for all values of x where v n is super-fast-magnetosonic. In the reconnection 
downflow, the fast-magnetosonic speed is effectively equal to the sound speed and we can write the 
shock formation condition as v n > Cf. 
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Fig. 9. — Overlaid plots of the dual constraints of the shock reprocessing model for pulse strengths 
and the shock formation model. Left S s = 4, Right: 5 S = 5. The inner box represents the physically 
allowed region. 
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Comparison to Aschwarideri Fits 
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Fig. 10. — r vs. E for the cooled and uncooled STFA models as well as the empirical E 3 / 2 fit of 
Aschwanden, et.al. (1995) 



